Collective dynamics in molten potassium: an Inelastic X-ray Scattering study. 
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The high frequency collective dynamics of molten potassium has been investigated by inelastic 
x-ray scattering, disclosing an energy /momentum transfer region unreachable by previous neutron 
scattering experiments (INS). We find that a two-step relaxation scenario, similar to that found 
in other liquid metals, applies to liquid potassium. In particular, we show how the sound velocity 
determined by INS experiments, exceeding the hydrodynamic value by ~ 30%, is the higher limit 
of a speed up, located in the momentum region 1 < Q < 3 nm~^, which marks the departure from 
the isothermal value. We point out how this phenomenology is the consequence of a microscopic 
relaxation process that, in turn, can be traced back to the presence of "instantaneous" disorder, 
rather than to the crossover from a liquid to solid-like response. 

PACS numbers: 67.55.Jd; 67.40.Fd; 61.10.Eq; 63.50.-(-x 



I. INTRODUCTION 

Alkali metals have been known since long time to ex- 
hibit remarkably pronounced inelastic features (Brillouin 
peaks) in their density fluctuations frequency spectra, 
even for wavevector (Q) well beyond the region of valid- 
ity of simple hydrodynamics. For this reason, in the last 
thirty years, the microscopic dynamics of such systems 
has been widely investigated both experimentally, by 
means of Inelastic Neutron Scattering (INS))^i^i^i'^i^i'^i'^i^ 
and numerically^i-^^i^^i^^i^'^i-^'^i^^i^^ with the aim of clar- 
ifying in detail the mechanisms underlying the atomic 
motions at a microscopic level. 

More recently, the possibility of studying experimen- 
tally the microscopic dynamics has been significantly in- 
creased by the advent of Inelastic X-ray Scattering (IXS) 
spectroscopj*i2*i^. Thanks to this technique, new insight 
into the collective properties of alkali metals have been 
provided'^^i^'-'i^^i^^i^^'^^, allowing to draw a much clearer 
picture of the high frequency dynamics, so that a deep 
comprehension of the ultimate nature of the high fre- 
quency excitations in these systems is gradually emerg- 
ing. In particular, thanks to development of approaches 
based on the memory function formalismfSS*^ a frame- 
work has been developed which accounts for the re- 
laxation spectra of the density fluctuations in a large 
wavevector regionSi. This framework applies also to 
non-alkali metalsSSiSS, hydrogen bonded liquids^Si^ and 
polymers^^. More specifically, it has been established 
that the decay of the density-density correlation function 
occurs through different relaxation mechanisms charac- 
terized by different time scales^i. These mechanisms ap- 
pear as additive contributions to the second order mem- 
ory function {AI{Q,uj)) of the density fluctuation spec- 
trum. Indeed, in addition to the thermal relaxation pro- 
cess, predicted by the hydrodynamic equations of a sim- 
ple fluid and arising from the coupling between density 
and temperature fluctuations, one must also allow ex- 



plicitly for a non Markovian behavior of the stress tensor 
correlation function thus introducing viscous relaxation. 
The most palpable consequence of these relaxation pro- 
cesses, is that the sound propagation is no longer deter- 
mined by the hydrodynamic sound velocity. 

On the basis of spectroscopic experimental evidence, 
viscous relaxation is thought to proceed throughout two 
distinct channels which are active on quite different time 
scalesSi. A first rapid decay of M{Q,uj) (on the time- 
scale Tfj^iQ)) is generally described as arising from the 
interactions of each atom with its "cage" of surround- 
ing neighbors. This is followed by a slower decay process 
(time-scale Ta{Q)) with long lasting tails which originates 
from slow, temperature-dependent structural rearrange- 
ment. Eventually, in those systems capable of under- 
cooling, these rearrangement undergo a structural arrest 
thus giving rise to the glass transition phenomenon. Ob- 
viously, in order to fully understand the characteristics of 
these processes, which are in principle common to many 
liquid systems, thorough investigations in an extended 
Q — E region are required. Indeed, in the recent past, 
massive investigations by means of inelastic neutron scat- 
tering (INS), devoted to liquid metals, have produced sig- 
nificant advances in the comprehension of the collective 
properties of liquids, in particular as far as single parti- 
cle properties are concerned. Nevertheless, with the aim 
of studying collective properties, these approaches have 
often suffered from limitations, intrinsic to INS, which 
arise mainly from the presence of incoherent scattering 
and from kinematic constraints restricting the accessible 
Q — E region. Very recently, accurate INS experiments 
have been performed^ on liquid potassium in which the 
inelastic components of the spectra have been analyzed 
either within a memory function framework*" or by em- 
pirical fitting functions^. In both these studies, values of 
sound speed exceeding the hydrodynamic one have been 
reported, and this has been interpreted as a reminiscence 
of crystal-like sound propagation. This interpretation is 
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substantially different from the results reported in other 
liquid metals (Li, Na, Al, Ga)S2i2ii22i2S and thus deserves 
a deeper investigation. 



II. THE EXPERIMENT 

In this work, we present an IXS study of the dynamic 
structure factor S{Q,uj) in liquid potassium. Thanks to 
i) the purely coherent nature of the IXS cross section 
in monatomic systems, which gives direct access to the 
collective dynamic structure factor, and ii) the lack of 
kinematic limitations, we have been able to derive con- 
vincing information on the ultimate nature of the atomic 
motion in this system. In particular, we show that, as 
reported in several other simple liquids^^ for Q values 
ranging from just below the diffraction peak maximum 
{Qm « 16 nm^^) all the way down to the hydrodynamic 
region (Q « 1 nm~^), the collective dynamics proceeds 
basically through two distinct viscous relaxation pro- 
cesses. Within this framework, the observed speed up of 
the sound velocity, reported by the INS measurements^i^, 
can be interpreted as due to a relaxation mechanism. In 
particular, thanks to the accessibility of the low-Q and 
high-i? regions, we show that the process mainly respon- 
sible for this increase in the sound velocity is not related 
to a crossover from a liquid-like to solid-like response. In- 
deed, although the system is above the melting point, the 
diffusive motion is frozen on the timescale of the density 
fluctuation even at the lower investigated Q values (i. e., 
in terms of the Brillouin peak frequency lobi one finds 
'-^ B{Q)Ta{Q) » 1 over the whole investigated Q range). 
On the contrary, the characteristics of the sound propa- 
gation are controlled by the faster relaxation induced by 
the peculiarity of the atomic vibrations associated with 
the instantaneous positional disordei^SiS. 

The experiment has been carried out at the ID16 
beam-line of the ESRF in fixed exchanged wave-vector 
configurations. A typical energy scan (—50 < E < 50 
meV) took about 300 minutes and was repeated several 
times in order to achieve a total integration time of about 
300 seconds per point. A five analyzer bench allowed us 
to collect simultaneously spectra at five different values 
of the exchanged wave- vector Q for each single scan. The 
sample thickness, (« 1 mm) was chosen in order to match 
the absorption length at the incident energy of 21748 eV, 
corresponding to the (11 11 11) reflection from the sil- 
icon analyzers. The energy resolution, 5E, in the used 
configuration was ~ 1.5 meV FWHM. The sample was 
kept under vacuum at a temperature of 343±1 K with- 
out need of any further container windows on either the 
beam incidence and viewing sides. 

In Fig. n we report the measured IXS intensity for 
each investigated (fixed) Q-value where the presence of 
an acoustic propagating mode clearly appears directly 
from the raw data. Also visible from the data is the pres- 
ence of a maximum of the dispersion relation (at Q 
nm~^) and the DeGennes narrowing at Q-values close to 
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FIG. 1: IXS spectra of liquid potassium (r=343 K) at the 
indicated Q values (circles) 



Qm, the position of the main peak of the static structure 
factor S{Q) (Q™ -16 nm-^). 



III. DATA ANALYSIS 

In order to extract quantitative information on these 
excitations, we performed a data analysis according 
to the generalized hydrodynami(j2Si2&, already tested 
against experimental data available on similar systemsSi. 
In this framework, the classical dynamic structure factor 
can be expressed in terms of its second order memory 
function, M{Q,t), so that: 



S{Q,u;) 



SiQ)n~^u;UQ)M'iQ,Lu) 



oj^ ^lj^{Q)+ujM"{Q,uj) + ujM'{Q,uj) 



where we have introduced the real {M'{Q,w)) and 
imaginary (M" (Q,w)) part of the Fourier transform of 
M{Q,t), and the frequency w^(Q) = KbTQ'^ /mS{Q) 
which is related to the generalized isothermal sound 
speed Ct{Q) — ijJvi{Q)/Q- This latter quantity can be cal- 
culated from the liquid structure once S{Q) is known (m 



3 



is the atomic mass and Kb the Boltzmann constant). In 
order to fit the data we have accounted for the detailed 
balance and for the convolution with the instrumental 
resolution fmiction R{uj) so that the quantity Im{Q,oj), 



Im{Q,^^) 



huj'/KT 



1 



-hui' /KT 



S{Q,uo')R{uo-uj')duj', 



is compared with the measured scattered intensity, 
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FIG. 2: Experimental S{Q,u)) (open circles) for selected Q 
values plotted together with the fitting function (full line) 
described in the text. The instrument resolution function 
(full- width at half hight SE « 1.5 meV) is also shown (dotted 
line) . 

Taking advantage of the results obtained from several 
other liquid metals, we utilized a memory function char- 
acterized by two relevant time-scales associated to two 
processes of viscous origin, besides the usual decay term 
accounting for thermal relaxation process. This scheme 
is able to account for the whole S{Q,uj) features, i.e. for 
both the quasielastic peak and for the Brillouin compo- 
nent. To represent the two viscous terms we adopted 
two simple exponential decays (corresponding to Debye 
relaxations) to account for the structural (a) relaxation 



and for the faster, microscopic (/i) contribution. The 
thermal relaxation has been described within the usual 
hydrodinamics result, which also predicts a Debye shape. 
Consequently, the total memory function reads 

M{Q,t) = (7-l)c.2(g)e-^-Q=* 

+ A2(Q)e-*/-"(Q)+A2(Q)e-*/v(Q). 

The value of ujo{Q) has been calculated using the S{Q) 
data reported in Refj^, while the specific heat ratio 7 and 
the thermal diffusion coefficient Dt are derived from lit- 
erature data"^"' neglecting their Q dependence. This lat- 
ter assumption is substantiated by the observation that 
in a very similar system, namely molten lithium, ^{Q) 
(and therefore the magnitude of the thermal relaxation) , 
varies less then 10% up to Qm^^ ■ Although in the large Q 
limit 7(Q) is expected to approach the limiting value of 
5/3, appropriate for a non-interacting monatomic sys- 
tem, any recourse to generalized-hydrodynamics argu- 
ments becomes increasingly doubtful for wavevectors dis- 
tinctly larger than Qm- 

All these parameters are kept fixed in the fitting proce- 
dure, while free fitting parameters are the structural and 
the microscopic relaxation times and r^, and the re- 
spective relaxation strengths, Aa{Q) and A^((5). Prior 
to the fitting, the experimental data have been put in 
absolute scale by using the first moment sum rule. 



= / S{Q,uj')da;' = S{Q), 



(1) 
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The tilde indicates the dynamical structure factor cor- 
rected by the detailed balance prefactor. In terms of the 
first spectral moments of the measured intensity (sub- 
script 1) and of the resolution function (subscript R) one 
can deduce that: 



^(Q)-^^o(i) /o(o) 



2m 



(1) 



Consequently, the normalized spectrum reads 



J I{Q,uj)duj 



A comparison between the best fitting line-shape and 
the experimental spectra is reported in Fig.[21for selected 
Q- values. 

It is interesting to test the applicability of the simple 
viscoelastic model (single relaxation) to the experimental 
data. As already noticed in several other systems, this 
oversimplified approximation does not account for the 
lineshape details at the level of accuracy reached by IXS 
(see Fig. ig. 



FIG. 3: Detail of IXS spectra of liquid potassium at Q=6.0 
and Q = 10 nm~^ . Open circles: experimental data. Thicker 
line: fit obtained by a model with two viscous relaxation chan- 
nels. Thinner line: fit results according to the usual viscoelas- 
tic model. 



One of the most relevant informations that can be ob- 
tained from an inelastic scattering experiment is the gen- 
eralized sound velocity, defined as ci{Q) — aj;((5)/Q2S*2&. 
The quantity lji{Q), corresponds to the position of the 
maximum of the longitudinal current correlation spec- 
trum given by Jl{Q,uj) = {uj'^ / Q^) S {Q , oj) which we 
have calculated using the classical model S{Q,uj) ob- 
tained from our best fit procedure. In Fig.^ the maxima 
of the current correlation spectra (full dots) are com- 
pared with the result of recent INS studies (stars and 
open circles As can be noticed, the IXS data signif- 
icantly extend at low Q the accessible kinematic region, 
allowing us to identify the transition from the hydrody- 
namic region to the high frequency, mesoscopic regime. 

It is worth to emphasize that the quantity which is 
related to the sound velocity is the maximum of the lon- 
gitudinal current correlation function. This is, for in- 
stance, the quantity which continuously evolves, at in- 
creasing wavevectors, towards the single particle regime, 
when distinct Brillouin modes are no longer visible in 
the dynamic structure factor. Moreover, in the frame- 
work of generalized hydrodynamics, it is the maximum 
of Jl{Q,i-^) which undergoes a transition between the 
well defined quantities uJo{Q) and ujociQ) describing the 
relaxed and unrelaxed mechanical regimes. 

In Fig. we report the apparent sound velocity 
(full dots). In the low Q region, the measured appar- 
ent sound velocity has a maximum {ci{Qmax) ~ 2240 
m/s at Qmax ~ 3 nm~^) where it clearly exceeds the 
isothermal value (~ 1790 m/s) deduced from ultrasonic 
measurements^i, a behavior (the so called positive dis- 
persion of the sound velocity) that is common with many 
other simple fluids. It is worth to recall here that, as al- 
ready pointed out in several IXS studies (see for example 
Refi^i), in liquid metals, owing to the very high thermal 
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FIG. 4: a) The maxima of the current correlation spectra 
Q,{Q) obtained from the spectra (full dots) compared with 
the result of recent INS studies (stars'" and open circles*). 
The dotted line is the hydrodynamic isothermal limit, b) 
Sound velocities deduced from the present study: appar- 
ent (•), co(Q) (dotted line) and Coo(Q) (full line) deduced 
by structural data. Also reported are the infinite frequency 
sound velocity Coc{Q) (A) and the unrelaxed sound velocity 
of the structural relaxation process Cooa(Q) (V) derived from 
the fit. 

conductivity, in the Q region probed here truly diffusive 
entropy mode are no longer possible. In other words, the 
linewidth (DtQ^) of the thermal diffusion mode exceeds 
the Brillouin frequency. As a consequence, the positive 
dispersion of the sound velocity consists of a speed up 
of the isothermal rather than the adiabatic value, as it 
happens, instead, in usual liquids. In the same figure we 
also report i) the generalized (Q-dependent) isothermal 
sound velocity uoq{Q)/Q (dotted hue), calculated using 
literature value for S'((5)24), that constitutes the low fre- 
quency limit of the sound speed, and ii) the infinite fre- 
quency sound velocity, Coo{Q) (full hne), related to the 
fourth moment sum rule of S{Q,uj)^ and numerically 
estimated using the Price-Singwi-Tosi pseudopotential — 
and a pair distribution evaluated by molecular dynamic 
simulation, according to the expressionSSiS^: 
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w m niQ^ J oz^ 

(2) 

The up triangles are the values of CodQ) deduced from 
the fit as coo(g) = VMQ)^ + A„(Q)2 + A^{Q)yQ. 
The down-triangles are the limiting (high frequency) ve- 
locity associated solely to the structural relaxation pro- 
cess: CooaiQ) = ^ojoiQ)"^ + Aa{Qy/Q, i. e. the velocity 
which is reached when the structural relaxation process 
is fully unrelaxed {ujB{Q)Ta{Q) » 1) or, in other words, 
the solid-like response. As it can be seen, the apparent 
sound velocity is always larger than CooaiQ), indicating 
the presence of a second relaxation process capable to 
drive the sound velocity from Cooa(Q) to Coa(Q). From 
this latter observation, and from the comparison between 
CooaiQ) ^■iid Coo(Q), it appears evident how the a pro- 
cess only plays a minor role in the full positive dispersion 
effect. On the other side, it is the faster, microscopic 
process (for which we find ojb{Q)'''^{Q) « 1 in the whole 
examined Q range), that pushes the velocity towards a 
fully unrelaxed regime (coo(Q)). 

It is worth to mention that, in both Refgi^ the pos- 
itive dispersion effect is fully ascribed to the transition 
between a liquid-like and a solid-like regime, namely to 
the structural relaxation process. This claim stems on 
the basis of the similarity of the sound velocity value 
of molten potassium with the value for the crystalline 
acoustic phonons along the [1 0] direction. This inter- 
pretation is in contrast with the results stemming out of 
the present work as well as with the analogous findings 
in lithium and sodiumSSiSi. In our opinion, as the pos- 
itive dispersion is a feature ubiquitous in liquid metals, 
but not in the crystalline counterpart, one should search 
his origin among the specific differences between the two 
phases: the diffusional motion and the instantaneous dis- 
order. As the first candidate is ruled out by the previous 
considerations (on the snap-shot timescale probed by IXS 
the liquid appears frozen), it seems natural to associate 
this effect to the positional disorder of the instantaneous 
configurations. This experimental fact is substantiated 
by the evidence of positive dispersion in the microscopic 
dynamics of simulated glassesiSiSi^. 

Coming back to Fig. the slight overestimate of 
Coo(Q) values as obtained by the fit (up-triangles) with 
respect to those derived from numerical evaluation has 
been already recognized as an effect of the exponential 
shape of M(Q, t) in the i ^ limi1jSL2i. 

A further interesting issue concerns the sound attenu- 
ation (Brillouin linewidth): once more the main actor is 
here the microscopic process, which crosses the resonance 
condition in the explored IXS Q- window (i^_b(Q)''7i(Q) al- 
ways remains of the order of one). On the contrary, the 

r 

"amount of viscosity" associated with the struc- 

tural relaxation does not contribute to the Brillouin 
width, as a consequence of the solid-like response all over 



the explored Q range (the Brillouin frequency always ex- 
ceeds the inverse structural relaxation time). 




Q (nm-^) 

FIG. 5; a) Values of the longitudinal viscosity rii{Q) as de- 
termined by the experimental data (•), dashed line is a guide 
to the eye. The hydrodynamic value is also reported (< — ). 
Inset: partial contributions to the rii{Q) due to the a (o)and 
(a) processes, respectively. Data for the structural process 
are reported for Q > 3 nm"'^ only, as for lower values the de- 
termination of Ta is unreliable due to finite resolution effects. 
The dashed line shows the compatibility of the /x process with 
a behavior. 

Within the hydrodynamic framework, the longitudinal 
viscosity can be generalized accounting for its Q— depen- 
dence and estimated as the area under the total memory 
function. In fact, the whole longitudinal viscosity stems 
from the sum of the structural and microscopic contribu- 
tions: 

TjiiQ)=piAlT^ + AlT,)/Q'. (3) 

We show in Fig. as the low Q limit obtained by ex- 
trapolation of the fitted data is in agreement with the 
hydrodinamic value rii{Q ^ 0) = 1.11 cP—. In the inset 
of Fig. O we plot the single terms of the viscosity asso- 
ciated with the two different processes A^Tq, A^r^; we 
observe that: i) the two contributions are of the same or- 
der of magnitude in the investigated Q range, ii) the mi- 
croscopic process exhibits a Q-dependence which is com- 
patible with a law consistently with the data already 
reported in other liquids and glasses21. 

IV. CONCLUSION 

We presented an experimental study of the collective 
high frequency dynamics in liquid potassium at the melt- 
ing temperature. Evidence for collective acoustic modes 
has been found in a Q region extending beyond the hy- 
drodynamic regime up to one half of the structure factor 



6 



main peak. Far from being a mere fitting exercise, a gen- 
eralized liydrodynamic analysis allows for a quantitative 
determination of relevant dynamical/thermodynamical 
parameters such as the sound velocity, the structural re- 
laxation times and longitudinal viscosity. Although the 
presented approach is based on the arbitrary choice of 
the functional form of the memory kernel, it allows one 
to point out how many relaxation channels are involved 
in the evolution of the collective dynamics and, more im- 
portant, to determine how they affect the dynamical and 
transport properties. 

The advantages of using a memory function based 
framework, therefore, are not in a better agreement with 
the experimental data but rather in the physical mean- 
ing of the involved parameters. This approach, there- 
fore, allows to gather information on the physical pro- 
cesses underlying the relaxation dynamics, and offers the 
possibility to measure physical quantities and their fi- 
nite lengthscale generalization. Through this study, for 
example, we have been able to show how the Brillouin 
linewidth is not related to the whole longitudinal viscos- 
ity (as it happens in the usual light scattering regime in 
ordinary liquids) but rather to the "microscopic" part 
of the viscosity. In other words, the damping mechanism 
active in the THz region does not stem from atomic rear- 
rangements (the viscous flow is frozen out) but is due to 
the dephasing of the density fluctuations consequent to 
the non-plane wave character of the vibrational motion 
around the quasi-equilibrium positions. 

Along the same way, the positive dispersion effect can 
be regarded as a disordered-driven relaxation process 
with the sound velocity inflection occurring at the res- 
onance condition. On the basis of the reported findings, 
therefore, the interpretation of the positive dispersion 
in terms of a crystal reminiscent behavior, claimed in 
Refsi^ can be questioned. The viscous to elastic tran- 



sition, indeed, well observed in the present study and 
in the whole alkali metal series, occurs at frequencies 
below the IXS window and is quantitatively negligible 
compared to the total sound velocity increase observed. 
In this respect, it is worth emphasizing how, in order 
to gain a quantitatively accurate determination of the 
sound velocity, one has to look at the maxima of the 
whole coherent current spectra, rather than arbitrarily 
distinguishing between a solely inelastic and a quasielas- 
tic contribution^^'^^. 

In conclusion, this study substantiates previous find- 
ings for the collective dynamics in other liquid metals 
such as Li, Na, Al, Ga^^i^^i^^i^^. This is an important 
indication of how, despite quantitative differences, the 
high frequency dynamics in simple fluids exhibit univer- 
sal features which go beyond system dependent details. 
Since on the observed time-scale the structure of the liq- 
uid is frozen {u!B{Q)Ta{Q) >> 1), one can think of the 
high frequency dynamics as that of a system with well de- 
fined equilibrium position (a "glass"). At the same time, 
the solid-like response alone does not account for the re- 
ported spectral features, as the details of the dynamics 
(sound velocity, microscopic relaxation times, viscosity) 
are mainly determined by the second, fast, relaxation 
process which is associated with the vibrations of the in- 
stantaneous disordered structure, in agreement with nu- 
merical works performed on similar systemsiS*^. 
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